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The  sea  state  of  the  Beaufort  and  Chukchi  seas  is  controlled  by  the  wind  forcing  and  the  amount  of 
ice-free  water  available  to  generate  surface  waves.  Clear  trends  in  the  annual  duration  of  the  open  water 
season  and  in  the  extent  of  the  seasonal  sea  ice  minimum  suggest  that  the  sea  state  should  be  increasing, 
independent  of  changes  in  the  wind  forcing.  Wave  model  hindcasts  from  four  selected  years  spanning 
recent  conditions  are  consistent  with  this  expectation.  In  particular,  larger  waves  are  more  common  in 
years  with  less  summer  sea  ice  and/or  a  longer  open  water  season,  and  peak  wave  periods  are  generally 
longer.  The  increase  in  wave  energy  may  affect  both  the  coastal  zones  and  the  remaining  summer  ice 
pack,  as  well  as  delay  the  autumn  ice-edge  advance.  However,  trends  in  the  amount  of  wave  energy 
impinging  on  the  ice-edge  are  inconclusive,  and  the  associated  processes,  especially  in  the  autumn  period 
of  new  ice  formation,  have  yet  to  be  well-described  by  in  situ  observations.  There  is  an  implicit  trend 
and  evidence  for  increasing  wave  energy  along  the  coast  of  northern  Alaska,  and  this  coastal  signal  is 
corroborated  by  satellite  altimeter  estimates  of  wave  energy. 
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1.  Introduction 

The  extent  of  seasonal  sea  ice  in  the  Beaufort  and  Chukchi  Sea 
of  the  Arctic  Ocean  is  changing  (Jeffries  et  al.,  2013).  This  paper 
explores  the  timing  and  location  of  the  annual  ice  minimum  and 
transition  to  refreezing  conditions,  with  application  to  the  sea  state 
over  the  open  water  portion  of  the  domain.  The  sea  state  is  set  by 
the  wind  forcing,  the  open  water  fetch  distance  available  for  wave 
generation,  and  the  duration  of  time  over  which  the  waves  can  ac¬ 
cumulate  energy  from  the  wind.  The  wind  forcing  is  episodic,  and 
thus  best  interpreted  as  probabilities  for  events  (i.e.,  storms).  The 
open  water  distance,  by  contrast,  has  a  much  smoother  signal  that 
is  dominated  by  the  seasonal  retreat  and  advance  of  the  sea  ice.  It 
is  the  combination  of  these  signals  that  determines  the  sea  state 
of  the  Beaufort  and  Chukchi  seas. 

Trends  in  the  Arctic  sea  ice  have  been  examined  by  many  pre¬ 
vious  studies  (e.g.,  Wadhams,  1990;  Wadhams  and  Davis,  2000; 
Stroeve  et  al.,  2005,  2008;  Simmonds  and  Keay,  2009;  Kwok  and 
Untersteiner,  2011).  Meier  et  al.  (2013)  show  that  in  recent  decades 
the  Arctic  sea  ice  cover  has  thinned  and  become  more  seasonal, 
such  that  the  total  area  covered  is  nearly  30%  less  at  the  an¬ 
nual  minimum  than  the  corresponding  mean  from  1979  to  2000. 
Stammerjohn  et  al.  (2012)  show  that  the  duration  of  the  summer 
open  water  season  since  1979  has  become  much  longer  in  the 
Beaufort  and  Chukchi  seas  due  to  an  approximately  1.6  months 
earlier  ice-edge  retreat  in  spring,  followed  by  an  approximately 
1.4  month  later  ice-edge  advance  in  autumn.  Stammerjohn  et  al. 
(2012)  also  find  inter-annual  links  to  the  reduced  ice  extent  which 
are  attributed  to  heat  fluxes,  especially  increased  duration  of  sum¬ 
mer  solar  heating,  coupled  with  an  overall  thinner  ice  cover. 

Coincident  with  the  delay  in  the  timing  of  the  autumn  ice  ad¬ 
vance,  there  is  a  trend  towards  stronger  autumn  storms  in  recent 
years  (Serreze  et  al.,  1993,  2001;  Zhang  et  al.,  2004).  The  combina¬ 
tion  of  these  winds  and  increased  open  water  distances  is  expected 
to  create  high  sea  states  (Francis  et  al.,  2011;  Francis  and  Vavrus, 
2012;  Vermaire  et  al.,  2013;  Thomson  and  Rogers,  2014)  and  in¬ 
crease  air-sea  fluxes  of  heat  and  momentum,  particularly  in  the 
Beaufort  and  Chukchi  seas  (e.g.,  Simmonds  and  Keay,  2009).  Some 
studies  have  connected  reduced  ice  cover  with  specific  storm  ac¬ 
tivity,  such  as  in  August  2012  (Simmonds  and  Keay,  2012;  Zhang 
et  al.,  2013;  Parkinson  and  Comiso,  2013).  Of  these,  Parkinson  and 
Comiso  (2013)  conclude  that  the  storm  reduced  the  September  ice 
extent  minimum  by  an  additional  5  percent.  This  relatively  small 
effect  suggests  that  high  sea  states  may  be  the  result  of  diminish¬ 
ing  sea  ice,  but  that  high  sea  states  are  not  yet  the  leading  cause 
of  diminishing  sea  ice. 

However,  there  is  some  evidence  for  feedbacks  between  ocean 
surface  waves  and  the  loss  of  sea  ice  (e.g.,  Asplin  et  al.,  2012). 
There  are  also  feedbacks  between  waves  and  ice  formation, 
such  as  the  rapid  freezing  that  occurs  when  waves  cause  pan¬ 
cake  ice  to  develop  (Wadhams  et  al.,  1987;  Lange  et  al.,  1989). 
Waves  are  both  associated  with  the  formation  of  pancakes  and 
attenuated  by  the  pancakes,  such  that  large  areas  of  the  ocean 
can  freeze  quickly.  Although  this  process  is  typically  associated 
with  the  Antarctic  ice-edge  or  the  Eastern  Arctic,  it  is  possible 
that  this  process  will  become  important  in  the  Beaufort  and 
Chukchi  seas  of  the  Western  Arctic.  For  example,  this  process 
is  already  common  in  the  Sea  of  Okhotsk,  which  is  relatively 
sheltered. 

Here,  we  set  aside  the  many  interesting  questions  of  wave-ice 
interactions  (e.g.,  Squire  et  al.,  1995;  Squire,  2007)  and  focus  in¬ 
stead  on  the  large-scale  patterns  of  the  sea  state  in  the  Beaufort 
and  Chukchi  seas.  In  particular,  we  examine  emerging  trends  in 
the  probability  of  high  sea  states  in  the  Beaufort  and  Chukchi  seas. 
The  recent  work  of  Wang  et  al.  (2015)  indicate  the  wave  heights 
are  increasing  slightly  and  wave  periods  are  increasing  strongly 


as  a  result  of  reductions  in  ice  cover  (as  opposed  to  changes  in 
the  winds).  We  examine  these  trends  and  the  autumn  ice  ad¬ 
vance  stage  in  particular.  Section  2  describes  the  data  products 
and  model  hindcasts  used  for  the  analysis.  Section  3  presents  the 
results,  using  a  full  climatology  of  ice  products  and  a  sub-set  of 
wave  hindcasts.  Section  4  discusses  the  findings  and  corroborates 
the  coastal  signal  with  satellite  altimeter  estimates  of  wave  trends. 
Section  5  concludes. 

2.  Methods 

Analysis  of  ice  and  sea  state  trends  uses  satellite  products  and 
model  hindcasts  from  an  area-preserving  domain  shown  in  Fig.  1. 
The  domain  is  a  rectangle  which  is  constant  in  area  with  latitude, 
such  that  the  range  of  longitudes  included  must  expand  north¬ 
wards.  The  domain  is  selected  to  cover  the  full  extent  of  the  sea¬ 
sonal  variation  in  sea  ice  cover  from  the  middle  of  the  summer 
(1  August)  to  the  late  autumn  (31  October).  The  analysis  that  fol¬ 
lows  uses  this  rectangle  and  is  restricted  to  the  months  of  August, 
September,  and  October. 

2.1.  Sea  ice  satellite  products 

The  analysis  of  sea  ice  area  coverage  used  the  NASA  Goddard 
Space  Flight  Center  (GSFC)  Bootstrap  SMMR-SSM/I  Version  2  quasi¬ 
daily  time  series  (1979  to  2014)  of  sea  ice  concentration  from  the 
EOS  Distributed  Active  Archive  Center  (DAAC)  at  the  National  Snow 
and  Ice  Data  Center  (NSIDC,  University  of  Colorado  at  Boulder, 
http://nsidc.org).  The  day  of  autumn  ice  advance  and  spring  retreat 
is  identified  for  each  gridded  (25  by  25  km  pixel)  location  and 
for  each  sea  ice  year  that  begins/ends  during  the  mean  summer 
sea  ice  minimum  (from  mid-September  to  mid-September).  When 
identifying  day  of  ice-edge  advance  and  retreat,  an  annual  search 
window  is  defined  such  that  it  begins  and  ends  during  the  mean 
summer  sea  ice  extent  minimum  in  mid-September.  Within  this 
interval,  the  year  day  of  ice-edge  advance  is  identified  as  when  sea 
ice  concentration  first  exceeds  15%  (i.e.,  the  approximate  ice-edge) 
for  at  least  five  days.  See  Stammerjohn  et  al.  (2012)  and  Comiso 
(2000,  updated  2015,  2010)  for  further  details. 

Sea  ice  type  was  estimated  by  scatterometer,  following  Gohin 
and  Cavanie  (1994)  and  Girard-Ardhuin  and  Ezraty  (2012),  with  the 
goal  of  examining  trends  in  the  relative  amounts  of  first-year  ice 
versus  multi-year  ice.  The  sea  ice  type  results  are  similar  using  the 
Envisat  altimeter,  following  Tran  et  al.  (2009). 

2.2.  Wind  reanalysis  product 

The  wind  and  ice  product  used  for  wave  hindcasting  is  ERA- 
Interim,  which  is  a  global  reanalysis  of  recorded  climate  observa¬ 
tions  over  the  past  3.5  decades  (Dee,  2011).  The  spatial  resolution 
of  the  data  set  is  approximately  80  km  (T255  spectral)  with  60  ver¬ 
tical  levels  from  the  surface  up  to  0.1  hPa,  and  the  grid  employed 
is  0.75  deg  resolution.  ERA-Interim  is  produced  by  the  European 
Centre  for  Medium-Range  Weather  Forecasts  (ECMWF).  The  tem¬ 
poral  coverage  is  four  time  steps  per  day.  The  10-m  wind  product 
is  used  to  estimate  the  wind  input  to  the  wave  model,  following 
the  latest  source  term  formulation  given  in  Ardhuin  et  al.  (2010). 

2.3.  Wave  model  hindcast 

Wave  evolution,  and  thus  the  development  of  a  sea  state,  is 
modeled  by  the  Radiative  Transfer  Equation,  as  follows: 

dE 

T-  V  •  ( CgE )  =  Swjncj  —  Sbrk  +  Snl  —  S[ce,  (1 ) 

where  E(co,  6)  is  the  directional  wave  energy  spectrum  and  cg  is 
the  group  velocity  (Masson  and  LeBlond,  1989;  Young,  1999).  The 
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Fig.  1.  Region  of  analysis,  (a)  Map  of  bathymetry  and  the  area-preserving  rectangle  defining  the  domain.  Green  colors  show  land,  (b)  Projection  of  the  domain  in  latitude 
and  longitude. 


equation  describes  the  temporal  and  spatial  evolution  of  waves 
as  an  energy  budget  in  frequency  co  and  direction  6.  The  deep¬ 
water  source/sink  terms  are:  input  from  the  wind  Swind,  dissipation 
via  breaking  Sbrk,  nonlinear  interactions  between  wave  frequen¬ 
cies  Snl,  and  interactions  with  sea  ice  Sice.  This  is  the  basis  of  all 
contemporary,  i.e.,  third-generation,  wave  prediction  models.  Here, 
we  use  the  WAVEWATCH-III  model  of  the  US  National  Oceano¬ 
graphic  and  Atmospheric  Administration  (NOAA)  (Tolman,  1991, 
2009)  with  a  16  km  resolution  polar  stereographic  grid  (Rogers 
and  Campbell,  2009)  for  the  entire  Arctic.  The  wave  model  also 
imports  ice  concentration  fields  from  the  ERA-interim,  which  are 
used  to  estimate  the  effects  of  sea  ice  on  the  waves  using  the 
Tolman  (2003)  scheme.  Regions  with  concentration  less  than  25% 
and  greater  than  75%  are  treated  as  open  water  and  land  respec¬ 
tively.  Partial  blocking  is  applied  for  intermediate  ice  concentra¬ 
tions. 

The  wave  model  hindcasts  are  performed  for  the  minimum  ice 
months  (August,  September,  and  October)  for  whole  Arctic  dur¬ 
ing  the  years  spanning  1992  to  2014.  A  more  detailed  analysis  is 
conducted  for  the  years  2004,  2006,  2012,  2014.  These  four  years 
bracket  the  modern  ice  conditions,  and  include  2012  as  an  extreme 
within  the  ‘new  normal’. 

Analysis  of  the  wave  model  output  within  the  defined  Beaufort 
and  Chukchi  domain  applies  a  threshold  definition  of  ice  concen¬ 
trations  less  than  0.15  in  defining  “ice-free”  areas.  The  percentage 
of  the  domain  determined  to  be  “ice-free”  according  to  this  thresh¬ 
old  is  tracked  in  time  for  each  hindcast.  Subsequent  analyses  use 
time  series  of  spatial  averages  from  the  ice-free  grid  cells,  in  par¬ 
ticular:  total  wave  energy,  f fEdOdco ,  the  wave  period  at  the  peak 
of  energy  spectrum,  Tp,  and  the  wind  stress,  r.  Analyses  also  use 
histograms  of  the  significant  wave  heights  Hs  from  all  ice-free  grid 
cells  and  all  time  steps  (i.e.,  no  spatial  or  temporal  averaging),  with 
the  conventional  definition 

EpgHj  =  ff  EdOdco.  (2) 

Finally,  an  evaluation  of  the  large-scale  potential  of  wave-ice  inter¬ 
actions  uses  the  normal  component  of  wave  energy  flux  incident 
to  the  ice-edge,  given  by 

f  =  J  Ecg-  me ,  (3) 

where  n  is  the  local  unit  vector  normal  to  the  ice-edge.  The  result 
is  the  total  rate  at  which  wave  energy  leaves  the  open  water  and 
enters  the  sea  ice  (i.e.,  the  boundary  of  a  control  volume).  Fig.  2 
shows  an  example  of  the  model  hindcast  and  application  of  Eq.  (3). 


2.4.  Satellite  altimeter 

Additional  wave  products  used  are  from  satellite  altimeters: 
the  entire  Envisat  record  (Queffeulou  and  Croize-Fillon,  2012)  and 
CRYOSAT  altimetry  from  the  NOAA  Laboratory  for  Satellite  Altime¬ 
try.  The  altimeter  data  were  quality  controlled  and  calibrated  ac¬ 
cording  to  Zieger  et  al.  (2009). 

3.  Results 

3.1  Ice  cover  results 

Trends  in  timing  of  ice  advance  were  determined  from  the 
passive  microwave  record  over  the  period  1979-2014  using 
the  method  described  in  Stammerjohn  et  al.  (2012).  Over  this  span, 
the  timing  of  the  autumn  ice  advance  has  become  significantly 
later  throughout  the  Arctic.  Fig.  3  shows  a  map  of  the  rate  of 
change,  in  days  per  year,  for  the  date  of  the  ice-edge  advance.  The 
most  pronounced  change  has  been  in  the  Beaufort  and  Chukchi 
seas,  where  the  statistically  significant  trend  is  1.4  days  later  per 
year,  with  a  similar  trend  towards  earlier  open-water  in  the  spring. 
The  trend  is  particularly  strong  near  the  northern  coast  of  Alaska 
and  the  Chukchi  shelf,  where  recent  years  have  almost  an  addi¬ 
tional  3  months  of  open  water  from  the  spring  to  the  autumn  (rel¬ 
ative  to  previous  decades). 

The  inter-annual  variability  of  this  signal  is  shown  in  Fig.  4, 
which  uses  a  spatial  average  of  the  ice-advance  date  over  the  de¬ 
fined  Beaufort-Chukchi  domain.  The  ice  advance  date  is  simply  the 
day  of  the  year  that  the  ice  covered  portion  of  the  domain  begins 
to  increase.  The  linear  trend  is:  0.41  ±  0.07  days  per  year.  Note 
however,  that  the  trend  over  the  whole  domain  is  modest  com¬ 
pared  with  the  coastal  portion  of  the  domain  (where  the  average 
trend  is  1.2  ±  0.2  days  per  year.  Although  2012  was  the  minimum 
ice  extent  by  area,  2007  is  actually  the  latest  timing  for  autumn 
ice  advance  in  the  record. 

The  changes  in  timing  and  ice  area  are  likely  related  to  the 
loss  of  multiyear  ice.  Ice  type  for  the  years  1999  to  2009  (using 
QuikSCAT)  and  2008  to  2015  (using  ASCAT)  is  shown  in  Fig.  5.  As 
seen,  in  the  domain  with  which  we  are  concerned,  the  extent  of 
multi-year  ice  has  decreased,  with  the  most  dramatic  retreat  in 
the  period  from  2005  to  2009.  Simultaneously,  the  extent  of  the 
first  year  ice  features  an  upward  trend.  Similar  results  can  also  be 
found  in  Maslanik  et  al.  (2007),  2011).  Based  on  satellite  measure¬ 
ments,  these  authors  concluded  that  the  sea  ice  in  the  Arctic  is 
becoming  younger  and  thinner,  represented  by  the  extensive  loss 
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Fig.  2.  Example  WAVEWATCH  III  hindcast  showing  significant  wave  heights  (color  scale),  wave  directions  (white  arrows),  ice-edge  (magenta  curve),  Beaufort-Chukchi  domain 
(white  outline  box),  and  ice-normal  energy  flux  time  series  (lower  panel).  The  red  dot  in  the  lower  panel  corresponds  to  the  time  of  the  wave  height  map. 
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Fig.  3.  Average  rate  of  change,  in  days  per  year  (contours  and  colors),  of  the  timing  for  the  autumn  ice  advance  in  the  Arctic.  The  most  notable  delay  in  ice  advance  is  in 
the  Beaufort  and  Chukchi  seas  (north  of  Alaska).  Trends  greater  than  ±  0.5  days  per  year  are  significant  at  the  0.01  level,  with  standard  error  determined  using  the  effective 
degrees  of  freedom  present  in  the  regression  residuals. 
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Timing  of  Sea  Ice  Advance 


Fig.  4.  Spatial  average  for  the  date  in  the  autumn  when  sea  ice  begins  to  refreeze  and  advance  southwards,  by  year.  The  solid  black  line  is  the  average  over  the  entire 
Beaufort  and  Chukchi  domain.  The  gray  dashed  line  is  the  average  within  the  coastal  perimeter  of  the  domain.  The  trends  are  shown  as  thin  lines. 


2000  2005  2010  2015 
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Fig.  5.  Multi  year  (solid  line)  and  first  year  (dotted  line)  sea  ice  extent  estimate  in  the  Arctic  for  March  since  2000  using  satellite  scatterometers.  QuikSCAT  sensor  estimates 
are  in  blue,  ASCAT  results  are  in  red  (Ifremer/CERSAT). 


of  perennial  multi-year  ice.  Similarly,  the  long-term  reduction  in 
sea  ice  thickness  in  the  Arctic  was  clearly  identified  by  Kwok  and 
Rothrock  (2009)  using  a  combination  of  submarine-  and  satellite- 
derived  thickness  measurements. 

Both  the  spatial  view  of  the  overall  trend  (Fig.  3)  and  the  tem¬ 
poral  view  averaged  over  the  domain  (Fig.  4)  indicate  that  in  re¬ 
cent  years  the  Beaufort-Chukchi  domain  has  more  space  and  time 
with  open  water  in  the  autumn.  Coupled  with  the  known  pattern 
of  strong  winds  in  the  autumn,  the  logical  expectation  is  for  the 
sea  state  to  increase. 

3.2.  Sea  state  results 

The  relationship  between  the  changing  autumn  ice  advance  and 
the  sea  state  is  evaluated  using  wave  model  hindcasts  of  the  late 


summer  and  autumn  from  four  years  that  span  recent  ice  condi¬ 
tions.  The  2004,  2006,  and  2014  ice  conditions  are  used  as  “typ¬ 
ical”  years,  and  2012  is  used  as  an  extreme  year  (with  minimal 
ice  extent  and  delayed  ice  advance).  This  extreme  year  (2012)  had 
anomalously  high  air  and  sea  surface  temperatures  during  the  au¬ 
tumn  months,  and  this  probably  to  the  observed  delay  in  the  ice- 
edge  advance  relative  to  other  years. 

Fig.  6  shows  the  time  series  of  area-averaged  ice  and  sea  state 
quantities  from  these  hindcasts.  The  percent  of  ice  free  area  in  the 
domain  (panel  a)  is  a  relatively  smooth  quantity  in  time,  because 
of  area-averaging.  In  contrast,  the  sea  state  quantities  of  wave  en¬ 
ergy,  peak  period,  and  wind  stress  (panels  b,  c,  and  d,  respectively) 
have  high  variability,  because  the  sea  state  is  event-driven  and  the 
autumn  storms  often  encompass  much  of  the  domain  (such  that 
area-averaging  does  not  smooth  the  signal). 
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Fig.  6.  Time  series  of  spatial  averages  over  the  Beaufort  and  Chukchi  Sea  in  hindcasts  of  four  selected  years:  (a)  open  water  fraction,  (b)  wave  energy,  (c)  wave  peak  period, 
(d)  wind  stress. 


The  evolution  of  ice-free  area  for  the  four  hindcast  years  is  con¬ 
sistent  with  the  timing  of  autumn  ice  advance  (Fig.  4),  although  it 
is  interesting  to  note  that  2006  has  a  similar  ice-advance  to  2004 
and  2014,  despite  much  less  ice-free  area  in  the  late  summer.  The 
ice  free  area  and  the  delay  in  ice  advance  are  both  notably  larger 
for  2012  than  the  other  years.  This  means  more  time  and  space 
were  available  for  the  generation  of  waves,  given  a  set  of  wind 
forcing  conditions.  However,  the  time  series  of  wave  energy,  peak 
period,  and  wind  stress  are  not  noticeably  different  between  2012 
and  the  other  hindcast  years.  Indeed,  the  ‘Great  Arctic  Cyclone’ 
of  August  2012  is  hardly  evident  in  this  analysis.  All  years  show 
a  consistent  increase  in  winds  and  waves  into  the  autumn.  The 
largest  event  energy  is  actually  from  the  year  with  the  least  ice- 


free  area  (2006),  though  it  did  have  the  strongest  wind  event,  as 
described  below.  This  event  was  an  intense  storm  near  the  coast  of 
Alaska,  with  hindcast  26  m/s  maximum  winds  and  8  m  significant 
wave  height.  This  highlights  the  importance  of  wind  forcing  in  de¬ 
termining  the  sea  state,  even  with  large  variations  in  ice-free  area. 
Since  the  area-averaged  wind  is  not  noticeably  different  between 
the  different  years  (other  than  the  particular  storm  of  Oct  2006), 
it  is  not  surprising  that  the  area-averaged  waves  are  not  noticeably 
different. 

However,  the  event-driven  nature  of  the  sea  state  is  best  ex¬ 
amined  probabilistically.  Histograms  and  fitted  Weibull  probability 
distribution  functions  are  used  to  identify  differences,  and  this  is 
where  the  effect  of  a  low  summer  ice  extent  minimum  followed 
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Fig.  7.  Normalized  histograms  of  the  significant  wave  height  at  all  grid  cells  and  all  time  steps  for  each  of  the  hindcast  years.  Normalized  probability  distribution  functions 
for  significant  wave  height  at  all  grid  cells  for  each  of  the  hindcast  years. 
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Fig.  8.  Normalized  histograms  of  the  significant  wave  height  at  all  ice  free  grid  cells  and  all  time  steps  for  each  of  the  hindcast  years.  Normalized  probability  distribution 
functions  for  significant  wave  height  at  all  ice  free  grid  cells  for  each  of  the  hindcast  years. 


by  a  late  ice-edge  advance  in  autumn  in  2012  is  very  apparent. 
Using  the  whole  domain  and  all  time  steps  of  the  hindcasts  ad¬ 
dresses  probability  of  a  given  sea  state  anywhere  in  the  domain, 
with  an  explicit  dependence  on  ice  cover.  Restricting  the  analysis 
to  ice-free  grid  cells  addresses  the  probability  of  a  given  sea  state 
anywhere  there  is  open  water,  with  an  implicit  dependence  on  ice 
cover.  In  the  figures  that  follow,  results  from  both  the  whole  do¬ 
main  and  the  ice-free  portion  are  presented. 

Fig.  7  shows  normalized  histograms  of  significant  wave  heights 
and  fitted  probability  distribution  functions  for  each  year  using  all 
points  in  the  domain.  The  results  are  skewed  by  the  high  num¬ 
ber  of  points  with  sea  ice  cover  (and  thus  zero  or  negligible  wave 
heights).  The  2012  distribution  differs  from  the  other  years,  with  a 
higher  mean  ((Hs)  ~0.6  m  versus  (Hs)  ~0.3  m)  and  longer  tail.  For 
example,  the  2012  results  have  an  almost  10%  chance  of  2  m  waves 
at  any  grid  cell,  compared  with  a  1%  chance  of  this  wave  height  in 
the  other  years. 

Fig.  8  shows  normalized  histograms  of  significant  wave  heights 
and  fitted  probability  distribution  functions  for  each  year  using 
only  ice-free  points  in  the  domain.  The  ice-free  results  across  the 
different  years  are  more  similar  than  the  full  domain  results,  but 
2012  still  shows  the  largest  mean  and  highest  probability  of  larger 
waves  (except  in  the  very  tail  of  the  distributions,  where  limited 
sample  sizes  make  differences  statistically  insignificant). 

Fig.  9  shows  normalized  histograms  of  peak  wave  period 
and  fitted  probability  distribution  functions  using  only  ice-free 
points  in  the  domain.  Consistent  with  the  results  of  Wang  et  al. 
(2015)  and  the  expectations  of  wave  maturity  over  larger  distances, 
there  is  a  shift  to  longer  period  waves  for  2012.  More  striking, 


however,  is  the  distribution  for  2006,  which  is  the  year  with  much 
less  ice-free  area  but  similar  ice-advance  timing  to  2004  and  2014. 
The  average  2006  peak  wave  period  is  shorter  and  the  distribu¬ 
tion  of  peak  wave  periods  is  wider.  This  suggests  that  open  wa¬ 
ter  area  may  be  more  important  than  the  length  of  the  open  wa¬ 
ter  season  in  determining  sea  state,  since  the  area  difference  for 
a  year  like  2006  persists  throughout  the  whole  season  and  ap¬ 
plies  to  multiple  storm  events  (whereas  a  delay  in  ice  advance 
might  only  be  relevant  to  the  wave  evolution  of  a  single  storm). 
For  all  years,  the  wave  periods  are  still  short  (Tp  ~  6  s)  relative  to 
other  oceans,  indicating  that,  despite  the  emergence  of  swell  in  the 
Beaufort-Chukchi  domain  (e.g.,  Thomson  and  Rogers,  2014),  the 
sea  state  of  any  given  ice-free  location  in  the  domain  is  still  domi¬ 
nated  by  local  wind  waves. 

Returning  to  the  question  of  wind  forcing,  Fig.  10  shows  nor¬ 
malized  histograms  of  wind  speed  and  fitted  probability  distribu¬ 
tion  functions  using  only  ice-free  points  in  the  domain.  Although 
there  are  minor  difference  in  the  mean  wind  speeds,  the  storm 
winds  that  drive  high  sea  states  (>10  m/s)  are  not  significantly 
different.  This  is  consistent  with  Wang  et  al.  (2015),  who  find  that 
variations  in  wind  forcing  are  insufficient  to  explain  the  trends  in 
the  waves. 

To  examine  the  complete  signal,  wave  model  hindcasts  for  ev¬ 
ery  year  from  1992  to  2014  are  analyzed  following  the  same  fit¬ 
ted  Weibull  probability  distribution  function  analysis  used  for  the 
four  years  examined  in  detail.  Fig.  11  shows  the  Weibull  scale 
and  shape  parameters  for  significant  wave  height,  peak  period, 
and  wind  speed.  The  scale  is  used  as  a  proxy  for  the  mean  value 
and  the  shape  is  used  as  a  proxy  for  the  standard  deviation 
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Fig.  9.  Normalized  histograms  of  the  peak  wave  period  at  all  ice  free  grid  cells  and  all  time  steps  for  each  of  the  hindcast  years.  Normalized  probability  distribution  functions 
for  peak  wave  period  at  all  ice  free  grid  cells  for  each  of  the  hindcast  years. 
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Fig.  10.  Normalized  histograms  of  the  wind  speed  at  all  ice  free  grid  cells  and  all  time  steps  for  each  of  the  hindcast  years.  Normalized  probability  distribution  functions  for 
wind  speed  at  all  ice  free  grid  cells  for  each  of  the  hindcast  years. 


year 

Fig.  11.  Trends  in  the  Weibull  fit  parameters  for  significant  wave  height,  peak  period,  and  wind  speed  over  the  wave  hindcast  years.  Diamonds  are  the  scale  parameter  and 
the  vertical  bounded  lines  are  the  95%  confidence  intervals  of  the  shape  parameter  divided  by  a  factor  of  ten  (for  visual  simplicity).  The  black  dashed  lines  are  the  estimated 
trend  lines  of  the  scale  parameter.  The  significant  wave  height  scale  has  a  trend  of  0.01  m  increase  per  year  and  the  peak  period  scale  has  a  trend  of  0.04  s  increase  per 
year,  both  of  which  are  statistically  significant  at  95%  confidence.  The  wind  speed  scale  does  not  have  a  significant  trend. 
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Fig.  12.  Time  series  of  the  total  energy  flux  incident  (normal  component)  to  the  ice-edge  within  the  Beaufort  and  Chukchi  seas  for  the  hindcast  years. 


around  that  mean.  There  are  statistically  significant  trends  at  the 
95%  level  for  both  wave  height  and  peak  period,  but  not  for 
wind  speed.  The  peak  period  signal  is  particularly  important,  since 
most  wave-ice  interaction  studies  have  found  a  strong  dependence 
of  wave  attenuation  on  wave  period.  Following  Wadhams  et  al. 
(1988),  the  trends  in  Fig.  11  imply  an  increasing  penetration  scale 
for  waves  entering  the  sea  ice,  such  that  longer-period  waves  are 
expected  to  propagate  several  kilometers  into  the  ice  under  recent 
conditions. 

4.  Discussion 

It  is  logical  that  larger  ice-free  areas,  which  are  persisting 
longer  into  the  autumn,  will  result  in  higher  sea  states  occurring 
more  often  in  the  Beaufort  and  Chukchi  seas.  The  wave  hindcasts 
presented  here  support  this  prediction,  and  the  robustness  of  the 
result  lies  in  the  distinctness  of  the  mechanism:  all  that  is  required 
to  increase  the  probability  of  higher  sea  states  is  more  ice-free 
area,  and  secondly,  longer  ice-free  duration,  not  more  storms  or 
increased  wind  forcing.  A  compounding  mechanism  is  storm  du¬ 
ration:  if  storms  of  similar  magnitude  simply  persist  longer  over 
open  water,  the  resulting  waves  will  be  more  mature  and  carry 
more  energy  flux. 

The  impact  of  an  elevated  autumn  sea  state  on  the  overall  Arc¬ 
tic  system  is  difficult  to  determine  without  detailed  understand¬ 
ing  of  wave-ice  interactions,  coastal  impacts,  and  changes  to  fluxes 
across  the  air-sea-ice  boundary.  This  is  further  complicated  by  the 
event-driven  nature  of  the  processes.  A  simplistic  approach  to  the 
wave-ice  question  is  to  examine  the  total  wave  energy  flux  in¬ 
cident  on  the  ice  (Eq.  (3)).  This  is  distinct  from  the  question  of 
overall  wave  activity  (and  associated  air-sea  fluxes),  because  an  el¬ 
evated  sea  state  in  the  region  does  not  affect  the  ice  unless  the 
waves  reach  the  ice.  Paradoxically,  as  the  ice-free  regions  expand, 
there  is  more  room  for  localized  storms  that  are  far  from  the  ice 
and  may  not  directly  affect  the  ice. 

Fig.  12  shows  time  series  of  the  total  integrated  wave  energy 
flux  arriving  at  the  ice-edge.  Similar  to  the  energy  results  (Fig.  6), 
the  values  are  similar  across  the  years  and  generally  increase  later 
in  the  autumn.  This  suggests  that  waves  may  be  more  important 
as  a  mechanism  to  alter  ice  advance  (via  the  formation  of  pan¬ 
cakes,  etc)  in  the  autumn,  rather  than  as  a  mechanism  to  alter 
ice  retreat  (via  fracturing)  in  the  summer.  This  is,  of  course,  re¬ 
lated  to  the  increased  ice-free  area  for  wave  generation  in  the  au¬ 
tumn.  The  present  results  are  inconclusive  in  terms  of  trends  in 
wave  energy  flux  arriving  at  the  ice-edge.  Although  2012  had  more 
wave  activity  throughout  the  domain,  the  overall  rate  of  wave  en¬ 
ergy  arriving  at  the  ice-edge  was  similar  to  other  years.  Still,  the 


August  2012  storm  is  notable  and  waves  may  have  enhanced  the 
well-documented  effect  of  the  storm  on  the  rest  of  that  year  (e.g., 
Parkinson  and  Comiso,  2013).  Such  feedbacks  and  the  role  of  wave 
directionality  are  the  focus  of  publications,  such  as  Stopa  et  al. 
(2016). 

Given  that  wave  energy  flux  is  a  conserved  quantity,  with  only 
minimal  dissipation  occurring  as  waves  propagate  in  open  water 
(e.g.,  Ardhuin  et  al.,  2010),  the  increased  wave  energy  inside  the 
domain  during  the  2012  season  can  be  assumed  to  increase  the 
flux  along  the  other  boundary:  the  northern  coast  of  Alaska.  The 
satellite  altimeter  results  in  Fig.  13  corroborate  this  suggestion. 
Fig.  13  shows  a  statistically  significant  increase  in  wave  energy 
along  the  coast  from  2007  onward,  compared  with  no  significant 
trend  (and  an  apparent  slight  decrease)  in  the  wave  energy  along 
the  ice-edge.  The  satellite  altimeter  product  is  scalar  energy  only, 
and  thus  it  is  not  possible  to  calculate  wave  energy  flux  (Eq.  (3)) 
for  a  direct  comparison  and  reconciliation  with  the  wave  model 
hindcasts.  Moreover,  the  satellite  product  is  not  uniformly  sam¬ 
pled  and  is  poorly  suited  to  the  Weibull  distribution  fitting  that 
was  used  to  identify  trends  in  the  preceding  sections.  We  thus  rely 
on  the  model  hindcasts  for  overall  trends  in  the  wave  climate  and 
discount  the  non-significant  trend  in  the  altimeter  analysis. 

This  implication  for  increasing  wave  energy  along  the  coast 
is  significant,  given  the  highly  erodible  nature  of  this  coastline 
(Overeem  et  al.,  2011).  Furthermore,  this  would  suggest  that  winds 
are  preferentially  directed  off-ice.  If  so,  wind-wave  generation  in 
partial  ice  cover  may  become  more  important  in  the  future  Arctic, 
when  the  seasonal  marginal  ice  zone  is  expected  to  be  more  ex¬ 
pansive.  The  process  of  wind-wave  generation  in  partial  ice  cover 
is  likely  far  more  complex  than  present  models  suggest  (Li  et  al., 
2015;  Zippel  and  Thomson,  2016)  and  is  in  acute  need  of  improved 
understanding. 

5.  Conclusion 

The  autumn  storms  that  regularly  occur  in  the  Beaufort  and 
Chukchi  Seas  are  elevating  the  sea  state  now,  and  will  continue 
so  into  the  future,  simply  because  it  is  increasingly  likely  that  the 
storms  will  occur  over  larger  open  water  areas  that  persist  longer 
into  autumn.  It  is  yet  to  be  determined  if  the  higher  sea  states 
will  in  turn  feed  back  to  the  large-scale  evolution  of  the  sea  ice. 
The  increasing  sea  state  may  affect  not  only  the  ice  cover  devel¬ 
opment,  but  also  wave  forcing  in  the  coastal  zone.  Either  way,  the 
increasing  sea  states  may  alter  air-sea  fluxes  and  associated  ecosys¬ 
tem  processes.  It  is  possible  that  the  increasing  sea  state  may  play 
an  important  role  in  modulating  the  presumed  changes  in  air-sea 
fluxes  and  upper  ocean  properties  that  are  occurring,  and  in  turn 
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Fig.  13.  Yearly  results  from  satellite  altimetry  estimates  of  spatially  averaged  wave  energy  along  the  northern  coast  of  Alaska  (red),  along  the  ice-edge  (blue),  and  over  the 
entire  domain  (grey).  Dashed  lines  show  calculated  trends.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of 
this  article.) 


Ice  volume  in  Oct  over  the  Beaufort-Chukchi  domain 


Fig.  14.  Estimates  for  ice  volume  in  the  Beaufort-Chukchi  domain  for  1970-2100  in  the  months  of  August,  September  and  October  using  coupled  ice-ocean  model  following 
the  IPCC  AR4  climate  change  scenario  A1 B. 
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may  modulate  the  response  of  sea  ice  to  climate  change.  Finally, 
higher  sea  states  are  of  operational  importance  to  mariners  and 
seabed  drilling  operators  in  the  region,  for  whom  higher  sea  states 
can  increase  the  likelihood  of  dangerous  icing  conditions  on  ships 
and  structures. 

New  observational  data  has  just  been  collected  to  assess  many 
of  these  processes:  the  Office  of  Naval  Research  “Arctic  Sea  State 
and  Boundary  Layer  Physics”  program  (Thomson  et  al.,  2013)  fol¬ 
lowed  the  ice-edge  advance  during  autumn  2015  while  simultane¬ 
ously  sampling  in  situ  air-sea-ice  interactions  from  the  R/V  Sikuliaq 
and  multiple  autonomous  platforms.  Pancake  ice  associated  with 
wave  forcing  was  ubiquitous  during  the  field  campaign,  and  the 
importance  of  this  ice  type  is  assumed  to  be  increasing  with  the 
wave  climate  in  the  region.  The  Sikuliaq  cruise  report  and  related 
information  are  available  at  http://www.apl.uw.edu/arcticseastate. 

Such  process  studies  are  essential  to  constrain  the  imperfect, 
yet  necessary,  parameterizations  used  in  climate  models.  Climate 
predictions  for  the  Beaufort-Chukchi  domain  already  indicate  that 
the  expansion  of  seasonal  open  water  will  only  accelerate  in  the 
coming  decades.  Fig.  14  shows  one  such  example  of  the  predicted 
dramatic  decrease  in  ice  volume  through  the  autumn,  using  cou¬ 
pled  ice-ocean  model  following  the  IPCC  AR4  climate  change  sce¬ 
nario  A1B  and  results  from  Long  and  Perrie  (2013),  2015).  These 
ice  predictions  are  consistent  with  AR5  results  following  the  recent 
work  of  Wang  and  Overland  (2015).  Incorporating  the  feedbacks 
associated  with  a  changing  sea  state  may  significantly  alter  these 
predictions,  but  that  remains  a  speculation  until  the  processes  can 
be  quantified  and  applied  within  the  climate  models. 
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